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We study the topologically non-trivial semi-metals by means of the 6-band Kane model. Existence 
of surface states is explicitly demonstrated by calculating the LDOS on the material surface. In 
the strain free condition, surface states are divided into two parts in the energy spectrum, one 
part is in the direct gap, the other part including the crossing point of surface state Dirac cone is 
submerged in the valence band. We also show how uni- axial strain induces an insulating band gap 
and raises the crossing point from the valence band into the band gap, making the system a true 
topological insulator. We predict existence of helical edge states and spin Hall effect in the thin 
film topological semi-metals, which could be tested with future experiment. Disorder is found to 
significantly enhance the spin Hall effect in the valence band of the thin films. 

PACS numbers: 



I. INTRODUCTION 

Topological insulators (TIs) have been recognized as 
novel states of quantum matter. 1 They are of tremendous 
interest to both fundamental condensed matter physics 
and potential applications in spintronics as well as quan- 
tum computing. TIs are insulating in the bulk which 
is usually due to direct bulk band gaps. However, on 
the boundary there are topologically protected gapless 
surface or edge states. Usually in these materials the 
spin orbital coupling is very strong such that the conduc- 
tion band and valence band are inverted. Such phenom- 
ena have already been reported in the literature of the 
1980'spHID However, the topological nature behind these 
phenomena was not revealed at the time. With the recent 
development of topological band theories, a series of ma- 
terials including HgTe quantum well, Bi^Sbi-^, Bi2Se3, 
and Bi2Te3 have been theoretically predicted and exper- 
imentally realized as 2D and 3D TIsPH^ The search for 
TIs has been extended from these alloys and binary com- 
pounds to ternary compounds. Very recently a large 
family of materials namely the Heusler-related and Li- 
based intermetallic ternary compounds have been pre- 
dicted to bepromising 3D TIs through first principle 
calculations. 16 22 The enormous variety in these com- 
pounds provides wide options for future material syn- 
thesizing of TIs. However, despite an insulating bulk is a 
critical prerequisite to the TI theory, none of these newly 
found materials are naturally band insulators. Many of 
them are semi-metals or even metals. It has been shown 
that the band structure and band topology of many of 
these compounds closely resemble that of the zinc blende 
structure binary compound HgTe and CdTeP^^ 

In this work we study the evolution of surface states 
in these topologically non-trivial compounds whose band 
structure closely resemble 3D HgTe by means of the 6- 
band Kane model. By studying the local density of states 
(LDOS) on the material boundary, we demonstrate ex- 
plicitly the existence of surface states in the direct band 
gap. Furthermore, we show that in the strain free con- 



dition, the surface states are separated into two parts in 
the band structure. One part exists in the direct band 
gap, the other part of surface states (including the cross- 
ing point of Dirac cone) submerge in the valence bands. 
The latter surface states have distinct momentum depen- 
dent spatial distribution from the former. By applying 
uniaxial strains, the crossing point can be raised up into 
the strain induced insulating gap. In the thin films made 
of these materials, topologically protected helical edge 
states will emerge on the sample edges. In this way we 
show a crossover between 3D TIs and 2D TIs. In the 
strain free condition, the bulk band gap can be controlled 
by tuning the film thickness. Meanwhile, we found that 
in the thin films disorders or impurities can significantly 
enhance the spin Hall effect (SHE) when the Fermi level 
is in the valence band. 
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FIG. 1: (color online) A schematic illustration of the band 
inversion between T6 (red curve) and T8(blue curve), the left 
is the normal case where the blue curve represents the LH and 
HH of Ts valence band, the right is the inverted case where 
the LH flips up and becomes the conduction band, the Yq 
appears below the HH band. 
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II. MODEL HAMILTONIAN 

3D HgTe and CdTe share the same zinc blende struc- 
ture. The band topology of these two materials is dis- 
tinguished by the band inversion at the T -point, which 
happens in HgTe but not CdTe. This causes HgTe t o be 
topologically nontrivial while CdTe is trivial .EE 11 ! 28 ! The 
essential electronic properties of both are solely deter- 
mined by the band structure near the Fermi surface at 
T-point, where the bands possess T 6 (s type, doubly de- 
generate), Ts (p type, j = 3/2 , quadruply degenerate) 
and (p type, j = 1/2, doubly degenerate) symmetry.^ 
The band inversion in HgTe takes place because the Tq 



bands appear below the band, whereas in the normal 
case (such as CdTe) T 6 is above T 8 (see Fig l)! 23|24 | In this 
work, we use a 6-band Kane model Hamiltonian which 
takes into account the Tq and Tg band. The spin-orbit 
split off T 7 band usually appears far below the T 6 and 
Tg bands and hence can be neglected because it doesn't 
affect the low energy approximation. The 6-band Kane 
model describes the band structure near the r-po int w ell 
and adequately captures the band inversion storyP^ As 
expected, we find gapless surface states exist in the direct 
gap of 3D HgTe while absent in CdTe. 

The 6-band Kane Hamiltonian writes: 
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TABLE I: 


Band structure parameters of HgTe at T 


= okP 




E P = 2m P 2 /h 2 F 7i 72 


73 


-0.3 eV 


18.8 eV 4.1 0.5 


1.3 



k-\- — k r- zlz ik- 



T = E g + B(2F + l)(kf l +k 2 z ), 
U = -B ll (kf + k 2 z ), 
V = -B-f{k 2 - 2k\\ 

where mo is the electron mass. For the simplicity of the 
physical picture, we have taken the axial approximation 
7 = (72 +73)/2 which makes the band structure isotropic 
in k Xl k y plane. E g , E pi F, 71, 72 and 73 are material 
specific parameters. The basis functions are denoted as 
(|^i), |^ 2 »r 6 , (|^ 3 >, hM, 1^5), |^6»r 8 . Here we take the 
parameters of HgTe just for illustration of the physics 
(Table. 1). The same physics should also happen in the 
recently discovered Heusler and Li-based ternary com- 
pounds which share similar band topologies^^. 



III. TWO TYPES OF SURFACE STATES IN 3D 
SEMI-METALS 

Without strain, HgTe is a semi-metal with an inverted 
band structure. At the V point, the conduction band 
(r 8 LH) touches the valence band (T$ HH) and the T 6 
is below the Ts band (as in Fig.l). In Hamiltonian Hq 
the band topology is solely determined by the parameter 
Eg. At the T point, LH and HH are degenerate at energy 
E = 0, and T 6 is at position E = E g . So when E g > the 
band structure is normal and insulating with a positive 
energy gap A = E g 



When Eg < the band structure 
is inverted Tq appears below LH and HH. In this case, 
the system is semi-metal since the bulk band gap A is al- 
ways because of the LH/HH degeneracy. We will show 
this semi-metal is topologically non-trivial by that gap- 
less surface states exists in the direct gap of LH and HH. 
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FIG. 2: (color online) (a) Surface LDOS of 3D HgTe without strain, bright line in the direct gap between LH and HH 
Fg bands indicates the first type surface state, bright regions in the valence band indicates the second type surface state, 
Etq < Eo < Ehh — Elh] (b) A insulating band gap is opened with strain T e = 0, U e = 0, V e = —0.224eV , the first and second 
type surface states become connected, Er 6 < Ehh = E < Elh', (c) same with (b) at T e = O.leV, U e = — 0.05eV, V e = — 0.25eV, 
note that the Tq band has been inverted with the HH Tg band and thus appears in the middle of Tg bands(LH and HH), i.e. 
Ehh < Etq < Eo < Elh- A highest LDOS limit is set and higher data points have been filtered to this limit for a clearer 
view of the whole spectrum. 



And by lifting the degeneracy of LH and HH by strain or 
finite confinement which makes A non-zero, the system 
instantly becomes a topological insulator. Whereas in 
trivial semi-metals, lifting the zero band gap only turns 
the systems into trivial band insulators. 

To study the topological properties of this Hamiltonian 
we transform it into a tight-binding model on a cubic 
lattice, where such approximation substitutions are used: 

ki —> -sin(kia), 
a 

*< ^(l-cos(k z a)), (2) 

here ki refers to k x ,k y and k z , a is the lattice constant 
which is taken to be 4 A in this work. This approximation 
is valid in the vicinity of T point. 

Surface states reside only on the system surface, which 
will project larger LDOS on the surface than the bulk 
states. Thus the surface LDOS can be studied to identify 
the existence of surface states. The surface LDOS is given 
by p(k) = —^TrlmGoo(k), where Goo is the retarded 
Green's function for the top layer of a 3D lattice. In a 
semi-infinite 3D system the surface Green's function can 
be obtained by means of the transfer matrix, 

^oo = (E — H 00 - ifoiT) -1 , (3) 
T = (E-H 00 -H 01 T)- 1 H^ (4) 

where Hqq and Hqi are the Hamiltonian elements within 
and between the layers (or super cells) and T is the trans- 
fer matrix. Usually EqQ can be calculated iteratively 
until T converges, which is quite time consuming. Here 
we use a fast converging algorithm proposed by San- 
cho et,el to calculate the transfer matrix™ In Fig. 2a we 
present the LDOS on an infinite xy surface, where z di- 
mension is semi-infinite. As expected, the LDOS clearly 



shows the existence of surface states between the LH and 
HH bands. Through further checking we find that these 
states indeed reside only on the surface boundary (see 
Fig. 3a). Its spatial distribution is in the decaying form 
beneath the surface. Interestingly, we also find another 
kind of surface state submerging in the valence bands 
(bright regions in the valence band of Fig. 3a) . It shows 
up between the inverted T 6 and HH T 8 bands. We con- 
firm its surface state nature by checking its spatial dis- 
tribution along the z direction (Fig. 3b and 3c). It is 
found that the spatial distribution of this surface state 
has very distinct form from that between the LH and HH 
bands. It bears the oscillating feature of LDOS of bulk 
states (see Fig. 3d) but is also clearly decaying beneath 
the surface. We call the former surface state as the first 
type and the latter the second type. Both types of sur- 
face state project much larger LDOS on the surface than 
the bulk states. However, the first type surface state is 
clearly decoupled from the bulk states while the second 
type is coupled with the valence bulk conduction states. 

The spatial distribution can be obtained by enlarging 
the super cell when calculating the surface Green's func- 
tion Gqo- For example, take the first 50 layers as a unit 
cell, we can obtain the LDOS distribution in the first 50 
layers. Fig. 3a shows the spatial distribution of the first 
type surface state at variant k. It shows that the closer 
to the T point the wider it distributes in space. Only 
away from the T point, it showes strong localization on 
the surface. Close to the T point, its distribution is bulk 
like and can barely be recognized as surface state. In 
Fig. 2 the bright line indicating first type surface has 
been highlighted for a clearer view. The Fermi surface 
of 3D HgTe is close to the poi nt wh ere the conduction 
band touches the valence bandP^ID The first type sur- 
face states at the Fermi surface distribute widely in the 
space. They are expected to make no significant contri- 
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bution to the transport properties either. Without ap- 
propriate doping or gating it is difficult to detect the 
first type surface states in 3D HgTe experimentally ei- 
ther through angle-resolved photo-emission spectroscopy 
(ARPES) or transport measurements. 

Unlike the first type surface states which show up as 
well defined sharp lines in the energy spectrum, the sec- 
ond type surface states appear as two crossing bright 
regions submerging in the valence bands in the energy 
spectrum. It appears near the boundary of the Tq band. 
At the crossing point, it is the sharpest, above the cross- 
ing point, it is broadened into two wide tails. Close to the 
T point, its LDOS on the surface is the highest. Away 
from the T point, the surface LDOS becomes smaller and 
finally those states merge into the bulk states. Also un- 
like the first type surface state whose distribution width 
becomes very large when approaching the T point, its dis- 
tribution width in space does not significantly depend on 
the momentum. At the T point, the second type surface 
state energy and wave function is exactly solvable from 
the Hamiltonian. We will show that this point is the 
crossing point of the surface state Dirac cone, (see text 



below) Beneath the surface, the LDOS quickly decays in 
an oscillatory way. For comparison, we also plot out the 
bulk states LDOS which is oscillatory and extended all 
over the space without decaying (Fig. 3d). Since from 
the LDOS the second type surface state can be easily dis- 
tinguished from the bulk states, we believe they may be 
easily detected from the ARPES measurements as well. 



IV. ORIGIN OF SURFACE STATES AND 
STRAIN INDUCED BAND GAP 

The origin of the two types of surface state can 
be understood intuitively in the following way. If we 
divide the Hamiltonian into two subspaces which are 
spanned by two groups of basis V\ = {^1, ^4, ^2, ^5} 
and V2 = {^3,^6}, it is amazing that the effec- 
tive Hamiltonian in the subspace V\ is quite similar 
to the effective model proposed for a 3D topological 
insulator.^ The Hamiltonian after arranging the basis 
as (^1(^1,^4,^2,^5)^2(^3,^6)) reads 
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Without considering the coupling between V\ and V2 
sub-spaces, V\ gives the LH (conduction) and Tq (va- 
lence) band, plus a gapless Dirac cone of surface states 
between their gap; V2 gives the HH (valence) which over- 
laps with the surface state Dirac cone and T 6 band. Af- 
ter turning on the coupling between the two sub-spaces, 
the surface state Dirac cone becomes separated into two 
parts. One part appears in the direct gap between LH 
and HH, the other submerges into the HH of valence 
band. In this sense, our results are very similar to an ear- 
lier work which discusses interface states in Hgi_ x Cd x Te 
hetero j unctions P 

The T 8 degeneracy at the Y point makes 3D HgTe a 
semi-metal. To lift the degeneracy and open an insulat- 
ing gap at the Fermi energy, which makes the system a 
true 3D TI, we consider applying a uniaxial strain along 



the (001) axis.E^ According to RefP, the additional 
strain induced Hamiltonian is introduced as (in the orig- 
inal basis order as in Eq. (1)) 








V 




T e 










u e + v e 











u e -v e 











u e -v e 



\ 



u e + v e ) 



where the T e , U e and V e are strain induced interaction 
terms. To illustrate the physical picture more clearly, 
we have made these terms artificially large in this work. 
At the T point, the boundary for the Tq band is Er 6 = 
Eg + T e , for HH band is E HH = U e + V e and for LH 
band is Elh = U e — V e . In the strain free condition, 
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FIG. 3: (color online) Surface state(ss) distribution along 
z direction. (a)First type ss at k=0.005(red), 0.01 (green), 
0.02(blue), 0.035(pink) 1/A; (b) second type ss below the 
crossing point at k=0.0066(red), 0.01 (green), 0.014(blue), 
0.017(pink), 0.022 (cyan) 1/A; (c) second type ss above the 
crossing point at k same as in (b); (d) bulk states at k=0.0066 
1/A and E=-0.4eV(red), -O.leV(green), 0.2eV(blue). 



Er 6 < Elh = Ehh — 0. However, under strain, the 
relative position of the three is subject to change and is 
determined by the strain terms. 

Surface state at the crossing point of the Dirac cone is 
easily solvable directly from the Hamiltonian Ho + H S1 
where k x = k y = 0. 30 Using the boundary condition 



ip(z = 0) = and ip(z = 
energy for surface state 



-oo) = 0, we get the eigen 



E 



where 



C 



(E g + T e )( 7 i + 27) + (U e - V e )(2F + 1) 
71 + 27 + 2F + 1 



(6) 



C 
M 



E g +T e + U e -V e 



E n 



U e + V e 



B 



(2F + 1-7X-27), 



B 



(2F + I + 71 +27). 



(?) 



And the wave function is 




with the basis ipi and ^4 or ^ 2 and ^5, where 



D 2 ± y/Dl - D 3 
2D+D_ 



with 



D 2 = -[Al + D+iEo-^ + D-iEo-I*)], 
D 3 = AD+D-iEo-L^iEo-^), 

D± = D x ± B u Li = C + M, L 2 = C - M and 



P. For the surface state solution to exist, it 



is required that: 



MB 1 > 0. 



This condition is easily satisfied in an inverted band 
structure without strain, where E g < 0. 

In Fig. 2b and 2c the spectrum after adding the strain 
interaction is plotted. An insulating gap is opened be- 
tween the LH and HH bands. Meanwhile, both the 
first and second type surface state go through an evo- 
lution. At a critical point (Fig. 2b) the crossing point 
can move to the top to the valence band at the V point 
and both types of surface states become connected with 
each other. On further increasing the strain strength, the 
crossing point jumps out of the valence band and sits in 
the band gap. At this stage, we find only one kind of 
surface state appearing in the insulating gap (Fig. 2c). 
The spectrum of surface state forms a gapless Dirac cone 
at the T point and the system becomes a typical 3D TI. 
Notice that Fig. 2c closely rese mbles the band structure 
of 3D TI Bi 2 Se3 and Bi 2 Te 3 ^^ Although Fig. 2 is 
obtained from the tight binding model on a lattice, we 
find the energy of the crossing point agrees with the exact 
solution quite precisely. 



V. EDGE STATES IN 2D THIN FILMS 
A. Quasi-2D lattice model 

It is however technically difficult applying strong 
enough strains to make a semi-metal insulating in ex- 
perimental conditions. The finite size effect can open a 
gap in the bulk states at the T point but will also open a 
gap in the surface states Instead of considering the 
surface states, we make the 3D topologically non-trivial 
semi-metals into thin films and study their edge effects 
in the strain free condition. The same lattice model is 
used as in the 3D case, but now z dimension is of finite 
thickness L. We study the LDOS on the side surface of 
a thin film where x is infinite and y is semi-infinite. In 
this case, only k x is good quantum number. In Fig. 4 the 
LDOS on the film edge surface for various film thickness 
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FIG. 4: (color online) LDOS at the edge of thin films at dif- 
ferent thickness calculated with 3D lattice model, (a) L=116 
A; (b) L=76 A; (c) L=36 A; (d) L=28 A; (e) L=20 A; (f) 
L=12 A. 



is plotted. In the thick film limit, the band gap opened by 
finite size effect is not obvious. The spectrum still resem- 
bles that of 3D case(Fig.4a and 4b). Surface states still 
show residuals on the spectrum. The finite confinement 
of z causes energy level discretization as in a quantum 
well, which induces a series of sub-bands appearing in 
the energy spectrum as layered structures. When the 
film is thinned down the discretized energy level spac- 
ing increases. A pair of edge states with linear disper- 
sion are found in the band gap(Fig.4c/d). The system 
becomes a 2D topological insulator. When the film is 
thinner than 20A, another transition happens, the sys- 
tem becomes trivial. The edge states disappear from the 
band gap. 



B. Quantum well approximation and 2D lattice 
model 

When the film is thin enough, the finite size caused 
band gap becomes obvious. In this case the finite confine- 
ment induced sub-bands are far away from the low energy 
regime. We can then use the quantum well approxima- 
tion (k z ) = 0, (k 2 ) ~ (7rAL) 2 p2 Using these relations in 
the Hamiltonian in Eq. (fiT) , and choosing the basis set in 
the sequence (|^i), I^Tfe), IV^), l^e), IV^)), we can 
obtain a two dimensional 6-band Kane model 

where 
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The system keeps time-reversal symmetry, and the rep- 
resentation of the symmetry operation in the new set 
of basis is given by T = /C • icr y (g) i3 X 3, where /C is 
the complex conjugation operator, a y and I denote the 
Pauli matrix and unitary matrix in the spin and or- 
bital space respectively. We can study the two blocks 
separately since they are time-reversal counterparts of 
each other. Here we focus on the upper block first. 
At k x = 0, the boundaries of T6,LH and HH are at 
E = E g + B(2F + l)(k 2 z ),E = -( 7l - 2-?)B(k 2 z ), and 
E = — (71 + 2j)B(k 2 ), which are controllable by choos- 
ing film thickness L. When L decreases from the thick 
limit, Down to L w 30^4 Tq band flips up and exchanges 



position with HH, the system is still non-trivial. Further 
down to L w 20A, Tq flips up and exchanges with the 
conduction band(LH)(see Fig.5e). The band structure 
then becomes trivial just as the illustrated case in Fig.l. 
This rough picture serves as an intuitive understanding 
of the topological transition and edge state formation in 
the thin films. It also agrees with result we obtained with 
3D lattice model in Fig. 4. In the thick film limit, the bulk 
band gap is always between LH and HH derived states 
while in thin limit it is between Tq and LH derived states. 
In the thick film limit, surface states (on top and bottom 
surfaces) and edge states(on side surfaces) co-exists. The 
first type surface states on the top/bottom surfaces be- 
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FIG. 5: (color online) LDOS at the edge of thin films at dif- 
ferent thickness using Hamiltonian (19k with 2D lattice model, 
(a) L=28 A; (b) L=25 A; (c) L=20 A; (d) L=16 A; (e) illus- 
tration of position change of LH, HH and T6 bands as film 
thickness is varied. 



come the effective bulk states of the film appearing as 
conduction band minimum. The edge states in the thin- 
ner films actually evolve out from the surface states on 
the side surfaces, whose nature is 2D instead of ID. 

Using the same approximation in eq. ^ , we can trans- 
form h(k) into a tight binding model on a 2D lattice. In 
Fig. 5 we show the LDOS on the edge a semi-infinite film 
for ft(k). The finite size gap agrees with that obtained 
with 3D lattice in Fig. 4 well. When L > 20A the edge 
states are found connecting the valence and conduction 
band. After the system becomes trivial when L < 20A, 
the edge states do not cross the band gap anymore, in- 
stead they only attach to the conduction band. At the 
critical point L = 20A, the valence band and conduction 
band touch and form a linear Dirac cone at the low en- 
ergy regime, which is showed in both Fig. 4 and Fig. 5. 
This shows that by controlling the film thickness, it is 
possible to obtain a single valley Dirac cone for each spin 
block without using the topological surface states .^1 No- 
tice that in Fig. 5 we can also see the edge states submerg- 
ing in the valence bands, which is similar to the case of 
second type surface states discussed previously. We call 
them the second type edge states. 



C. Chern number description of thin film band 
topology 

In Quantum Hall Effect (QHE), bulk-edge corre- 
spondence tells us that nonzero Thouless-Kohmoto- 
Nightingale-den Nijs (TKNN) integer, summed by the 
Chern number of occupied bands, is closely related to the 
presence of edge states on the sample boundaries! 33 1 34 1 For 
the time-reversal invariant systems which belong to the 
universality class of zero charge Chern number, Z2 index 
is then introduced to characterize the topologically non- 
trivial states!^ However, when the Hamiltonian of the 
system is composed of time-reversal diagonal blocks,^ Z2 
invariant can be identified with the parity of the Chern 
number for each block. 11 Focusing on the upper block of 
the thin film Hamiltonian Eq([8|(i.e h(k)), we can discuss 
its QHE. To identify the existence of edge states, we can 
calculate the first Chern number of each band and sum 
up for all occupie d ban ds. The Berry curvature for each 
band is defined a d 35 ' 36 1 



ft n (k) 



du n . 



dk x 



dk y 



dk v 



dk x 



Then the Chern number is computed by integrat- 
ing Berry curvature over the first BZ: C n = 
2^ J d 2 kQ n (k). For each block of the 2D lattice model, 
Chern number can be given to each of the three bands. 
Tuning the thickness of the quantum well L can inverse 
the band structure near the T point, resulting in the 
change of Chern number for the touching bands. Sev- 
eral critical thickness for L can be found when there's 
band touching at the T point: L w 30A and L « 20A. 
Numerical result of this integral divides the parameter 
space into three regions: for L G (30A, +00), Chern num- 
ber C = 1,0,-1 arranged in sequence from the lowest 
energy band to the highest (the same for the latter two 
cases); for L G (20 A, 30 A), Chern number C = 2, -1, -1 
and for L G (0, 20A), Chern number (7 = 2, -2, 0. When 
Fermi level lies between the highest and the middle en- 
ergy band, the TKNN integer TV = 1 for L G (20 A, +00) 
and N = for L G (0, 20 A). The former case is topolog- 
ically nontrivial( referred to in Fig. 5a and 5b) and the 
latter is trivial (referred to in Fig 5d) 



VI. QSHE AND DISORDER ENHANCED SPIN 
HALL EFFECT 

The edge state picture in thin films can beveri- 
fied with explicit Landauer-Biittiker calculations p3Ml We 
construct a four-terminal Hall bar device using the tight 
binding version of Hamiltonian (J8|. Because the sys- 
tem keeps time reversal symmetry (TRS), we need only 
consider a single block of the 6x6 Hamiltonian. The 
transverse charge conductance of each sub block is de- 



fined as G f 



G 



12 



Gin denotes the conduc- 



tance between terminal i and j, which is given by Gij 
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FIG. 6: (color online) Four terminal Landauer-Biittiker cal- 
culations for thin films. Four identical semi-infinite leads are 
attached to the cross bar device, which is shown schematically 
in the insert. Disorders are put on the central area(shaded ar- 
eas) . 200 disorder samples are taken for each disorder strength 
W. The arrow indicates the position of finite size gap of the 
helical edge states. 



^Tr[TiG R TjG A ]. r< = i[£f -£f] is the spectrum 
function of lead i. G R I A is the retarded and advanced 
Green's function of the sample which has taken into ac- 
count the four semi-infinite leads through the self energy 
E f/A : qr/a = [E -H c - Eti Ef ,A \- 1 - H c denotes 
the Hamiltonian of the shaded region of the sample.^®! 
In real system with TRS, the transverse charge conduc- 
tance calculated here corresponds to the transverse 'spin' 
conductance, which is defined as G ts — G\ c — G\ c . The 
longitudinal conductance is defined as G\ = G\ 3 + G\ 3 . 
However, it should be noted that the 'spin' up or down 
here is only an indice of each Hamiltonian sub block and 
does not mean real spin. Because the basises of each sub 
block are hybrid of j = ±1/2 and j = ±3/2 states of T 6 
and Ts bands P 

In Fig. 6, we show the result for a film whose thick- 
ness is 25 A. In the clean limit, i.e. without disorder, 
in the finite size induced band gap the G\ vanishes, Gt s 
is quantized as G ts — 2e 2 /h (with G\ 2 = G\ A = e 2 /h, 
G\ 4 = G\ 2 = 0) , which proves the existence of helical 
edge states in the thin films and indicates the existence 
of quantum spin Hall effect (Q SHE) ! 41 " 45 l In the conduc- 
tion band, the edge states coexist with the bulk states. 
G ts there is not quantized, which implies that the bulk 
states themselves carry transverse spin conductance that 
partly cancels out the quantized spin Hall conductance 
carried by the edge state channels. In the valence band, 
similarly, the second type edge states coexists with bulk 
states. Notice that the confinement in x,y dimension also 
generates a finite size gap in the edge state itself, which 
is indicated by the arrow in Fig. 6 where Gt s vanishes. 



Disorder is recently known to play an important role 
in some exotic phenomena in TIs! 46 * 47 ^ In real materi- 
als, dirty impurities are inevitable. It is even possible 
to control them artificially ! 48 * 49 ! Here we consider the ef- 
fect of disorder impurities to the transport properties of 
the thin films. Anderson type white noise is introduced 
on the lattice model, which are spin independent on site 
random potentials in the range [— W/2, W/2]. TRS is 
not voilated by disorder. In Fig. 6 we show the transverse 
spin conductance G ts of the four terminal device at vari- 
ous disorder strength W. The quantized Gts i n the band 
gap shows robustness against modest disorder strength, 
which is expected because the disorder we apply does 
not couple the TRS sub-blocks of the Hamiltonian and 
thus brings no backscattering between different spin edge 
channels!^ Interestingly, Gts is significantly enhanced 
in both conduction and valence band even with weak 
disorder. While G ts quickly becomes suppressed in the 
conduction band when disorder increases, it still remains 
enhanced in the valence band. Notice that with modest 
disorder strength(3eV in Fig. 6), the longitudinal conduc- 
tance G\ already vanishes in the low energy regime, but 
the transverse spin conductance remains enhanced in the 
valence band. Therefore with modest disorder strength, 
we expect a strong disorder enhanced spin Hall effect in 
the valence band, which would manifest itself through 
strong signals in non-local measurements. In a very re- 
cent experiment similar observation of SHE difference 
in valence and conduction band has been reported with 
HgTe quantum wells 



VII. SUMMARY 

We have demonstrated that 6-band Kane model can 
be utilized to study the topological properties of both 
3D and thin film realistic materials. By choosing proper 
model parameters, it describes the band structure well 
in the low energy regime and captures the essential 
band topological features sufficiently. By calculating the 
LDOS on the system boundary, the existence of topolog- 
ical surface states is explicitly demonstrated. Using the 
model parameters of HgTe, we show that even though the 
system is semi-metal surface states already exist in the 
material as long as the band structure is inverted. We 
also find in the strain free condition surface states are 
divided into two parts in the spectrum, each of different 
characteristics. We demonstrate that uniaxial strains can 
generate an insulating gap and transform the semi-metals 
into true TIs, in which the gapless Dirac cone of surface 
states is found in the bulk band gap. Because of the 
similar band structures and band topologies, we expect 
that the same physics applies to the recently discovered 
Heusler-related and Li-based intermetallic ternary com- 
pounds, which are mainly topologically nontrivial semi- 
metals or metals. 

We also demonstrate a crossover from 3D topological 
semi-metal to 2D QSHE insulators. In the thin films 
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made of these materials, we predict the existence of he- 
lical edge states and QSHE in the strain free condition. 
And finally we show that disorder plays an important 
role in the transport properties in the thin films. It sig- 
nificantly enhances the SHE in the valence bands. 
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